clc; %清楚前面的指令
close all; %关闭matlab所有窗口
clear all; %清空工作区
I=10;%初始患者
N=10000;%总人口
S=N-I;%易感人数
lemda=0.1;%接触率
mu=0.05;%治愈率
t=1:365;%时间
cgm_fu=mu/lemda;
for i=1:(size(t,2)-1)
    I(1+i)=I(i)+I(i)*(N-I(i))*lemda/N-mu*I(i);
    %I(1+i)=N/(1-cgm_fu) + (I(1)*(1-cgm_fu)*exp((lemda-mu)*i))/((1-cgm_fu)-I(1)/N);
    S(1+i)=N-I(1+i);
end
plot(t,I,t,S)
xlabel('时间')
ylabel('人数')
legend('患病者','易感者')
title('SIS传染病模型')